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Abstract 

An  study  of  the  one-dimensional  isentropie 
compression  experiment  (ICE),  performed  with  High 
Explosive  Pulsed  Power  (HEPP),  has  demonstrated  that 
accurate,  high  stress,  isentropie  Equations  of  State  (EOS) 
data  may  be  obtained  with  this  technique. 

The  physics  and  accuracy  of  electromagnetic  loading  in 
the  ICE  technique  is  presented.  It  is  shown  that  the 
HEPP-ICE  load  configuration  is  capable  of  producing 
magnetic  stresses  that  are  uniform  to  1  part  in  1000  over 
the  central  87%  of  the  sample  faces,  and  that  HEPP-ICE 
provides  exact  matching  of  the  stresses  between  opposing 
samples.  This  magnetic  uniformity,  the  exact  matching, 
and  the  large  sample  samples  possible  with  HEPP-ICE, 
are  necessary  for  the  highest  accuracy  isentropie  EOS 
data. 

Isentropie  EOS  data  have  been  obtained  with  a 
prototype  HEPP-ICE  system,  and  the  results  for  tungsten 
and  copper  demonstrate  the  inaccuracy  of  the  technique, 
which  may  be  as  low  as  0.2%  in  stress.  It  is  shown  that  a 
large-scale  HEPP-ICE  system  is  capable  of  producing 
shock-free  loading  up  to  2.2  TPa  in  10-mm  thick  tungsten 
samples. 

I,  INTRODUCTION 

The  basic  techniques  of  HEPP-ICE  have  been  described 
previously  [1,2,3].  The  original  ICE  technique  is  due  to 
Asay  [4],  and  a  significant  body  of  ICE  work  has  been 
conducted  on  the  Z-machine  at  the  Sandia  National 
Laboratory.  In  the  ICE  experiment,  smoothly  rising 
(shock-free)  mechanical  compression  waves  are 
propagated  into  matched  samples  of  different  thicknesses 
by  electromagnetic  loading  in  a  planar  geometry.  A 
complete  EOS  isentrope  is  acquired  in  one  experiment, 
i.e.,  continuously  from  zero  up  to  the  peak  stress.  Good 
quality  isentropie  EOS  data  have  been  obtained  for  many 
materials  [5,6,7]. 

A.  ICE  method 

Two  identical  samples,  with  a  difference  in  thickness  of 
h,  are  compressed  by  identical  B-forces  and  their  particle 


velocity  profiles,  u(t),  are  obtained  from  VISAR 
measurements  at  the  rear  free  surfaces.  Figure  1.  The 
magnetic  stress  is  given  by  the  vector  cross-product  Pb  = 
J  X  B,  where  J  is  the  current  per  unit  width  (dl/dW).  In 
the  simplest  case,  when  the  plate  width,  W  »  d,  the  plate 
separation,  the  magnetic  field  on  the  inside  surface  of  one 
conductor  due  the  current  in  the  opposing  surface  is  ‘A Bo  J 
(SI  units  are  used  throughout  this  paper).  Then  the  stress 
normal  to  the  inside  surfaces  Pb  ^  J  x  B  =  ‘AjuU^-  (As  d 
increases  with  time  and/or  when  the  current  is  non-steady, 
then  B  <  ‘AjUoJ^  and  becomes  non-uniform,  see  below.) 

Using  Lagrangian  wave  analysis  [8]  the  Lagrangian 
wave  speed  is  Cl(u)  =  At.  The  differential  form  of  the 
Lagrangian  momentum  conservation  equation  is 

da  =  pQCi(u)du  (1) 

which  is  used  to  calculate  the  change  in  stress,  da,  for 
each  change  (step)  in  particle  velocity,  du,  going  up  the 
curve  u(t),  where  po  is  the  initial  density  {h  and  po  are 
constant  in  Lagrangian  space).  Continuous  EOS 
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Figure  1.  Top:  Velocities  Ui  and  U2  are  measured  at 
the  outer  surfaces  of  samples  SI,  S2.  Bottom:  velocity 
Ci(u)  is  obtained  from  U ;  and  U2,  point  by  point. 
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relationships  between  stress,  wave  speed,  partiele 
veloeity,  ete.,  from  zero  to  the  peak  stress  are  thus 
obtained. 

B.  HEPP-ICE  apparatus 

A  6  mF,  20  kV  eapaeitor  bank  provides  a  seed  eurrent 
of  ~2  MA  to  a  4-inx5-in  plate  flux  eompressor  [2].  The 
FCG  transfers  up  to  12  MA  into  a  storage  induetor  of 
~25  nFl  and  an  explosively-formed  fuse  (EFF)  opening 
switeh.  Three  parallel  explosively-driven  polyimide 
elosing  switehes  [9]  then  transfer  eurrent  to  the  load  at  the 
appropriate  time.  This  prototype  eireuit  is  eapable  of 
delivering  7  MA  (dl/dt  ~3  xlO'^A/s)  into  loads  of 
typieally  1  to  2  em  width  and  will  produee  isentropie 
eompression  at  stresses  in  the  range  of  0  to  250  GPa. 
Flowever,  advaneed  HEPP-ICE  systems  are  eapable  of 
produeing  isentropie  data  beyond  2  TPa  [10]. 

C.  Stress  uniformity  and  exact  B-Jield  matching  in  an 
HEPP-ICE  load 

It  is  emeial  that  the  samples  are  subjeeted  to  identical 
and  uniform  magnetie  stresses.  Previously  we  showed 
that  for  statie  eurrent  flow  the  B-field  beeomes  inereasing 
non-uniform  aeross  the  width  as  the  ratio  of  separation  to 
width,  d/W,  is  inereased  [11].  However,  during  the  EOS- 
data  gathering  portions  of  ICE  experiments  the  eurrent  is 
always  rising  in  the  load,  i.e.,  dl/dt  >  0.  The  rise  is 
approximately  sinusoidal,  rising  from  0  to  7i/2  over  -'Aps. 
For  parallel  plates,  the  eurrent  flows  preferentially  along 
the  outside  edges  of  the  load,  where  the  speeifie  eleetrieal 
impedanee  is  a  minimum.  The  dynamie  eurrent  and  stress 
distributions  have  been  modeled  with  a  finite  element 
partial  differential  equation  eode  [12];  the  stress  plot  is 
shown  in  Figure  2. 


Figure  2.  Magnetie  stress  Pb  on  the  inside  surfaee  of  a 
parallel  eonduetor.  The  width  of  the  eonduetor  is  along 
the  A-axis,  the  thiekness  along  the  F-axis  (drawn  to  a  4x 
larger  seale  than  the  X  axis  for  elarity). 

The  ealeulation  shows  that  for  the  parallel  plate 
eonfiguration  used  in  HEPP-ICE,  under  typieal  dynamie 
eonditions,  the  stress  is  uniform  to  1  part  in  10^  aeross  the 
eentral  87%  of  the  eonduetors.  (The  plot  is  for  a  d/W  ratio 
of  0.3,  whieh  would  oeeur  towards  the  end  of  experiment, 
when  uniformity  is  at  its  worst.) 


Figure  3.  Left:  the  HEPP-ICE  load.  Right:  the  four 
samples,  SI  to  S4;  the  magnetie  fields  (@)  eommon  to 
samples  S1&S2,  B^;  and  the  eorresponding  field  for 
S3&S4,  B34. 

For  exaetly  matehed  magnetie  stress  loading,  the 
sample  positions  are  important.  In  Figure  3  two  opposing 
sample  pairs,  SI  and  S2,  S3  and  S4  are  shown  and  they 
share  the  magnetie  fields  B12  and  B34.  Currents  flow  in 
and  out  from  the  souree  at  the  bottom,  and  as  they  do  so, 
the  currents  redistribute  themselves  as  the  samples 
separate  under  magnetic  loading.  The  samples  will  have 
different  thicknesses  and  may  have  different  acoustic 
impedances  and  will  consequent  separate  with  different 
accelerations.  We  have  found  that  the  stresses 
corresponding  to  B12  and  B34  may  differ  by  as  much  as 
5%.  The  advantages  of  the  parallel  plate  configuration 
used  in  HEPP-ICE  are  that  the  B-fields  and  stresses  are 
exactly  matched  for  opposing  samples  because  they 
always  share  the  same  B-field  and  the  central  B-field  is 
very  uniform. 

D.  Higher  stresses,  larger  sample  sizes  and  higher 
accuracy  in  HEPP-ICE  experiments 
One  important  advantage  of  the  HEPP-ICE  over  other 
ICE  techniques  is  that  the  experiment  can  be  scaled  to  any 
magnitude  of  current,  and  hence  stress,  by  simply  by 
choosing  the  appropriate  HEPP  components.  When 
designing  an  ICE  experiment  it  is  important  to  prevent  the 
ramp  wave  developing  into  a  shock  wave  by  the  so-called 
shock-up  process.  Shock-up  is  caused  by  the  fact  that  the 
compression  wave  velocity  is  stress  dependent,  and  in 
most  materials  the  wave  velocity  increases  with  stress. 
Should  shock-up  occur  in  any  of  the  samples  in  an  HEPP- 
ICE  experiment,  the  compression  of  the  material  would 
no  longer  be  isentropie  after  the  shock  has  developed. 
The  maximum  stress  and  sample  size  ultimately  depend 
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on  how  well  the  shoek-up  in  a  material  ean  be  delayed. 
The  method  of  eharaeteristies  [13]  ean  be  used  to  design 
the  optimum  eurrent  profiles  to  prevent  shoek-up,  they  are 
shown  for  2.2-TPa  HEPP-ICE  experiments  and  20-mm 
wide  tungsten  samples  in  Figure  4.  Flere,  the  maximum 
thiekness  without  shoek-up  is  10  mm,  the  maximum 
eurrent  is  65  MA,  and  the  risetime  is  2  ps.  The  eurrent, 
peak  dl/dt  (40  TA/s)  and  risetime  are  well  within  the 
eapabilities  of  existing  HEPP  eomponents  presently  used 
at  LANE. 


within  1%,  whieh  is  also  the  estimated  error  in  these 
tungsten  experiments  eaused  by  the  diffieulties  of 
measuring  the  sample  thieknesses  to  the  required  aeeuraey 
in  this  partieular  experiment. 


Figure  5.  Tungsten:  pressure  vs.  partiele  veloeity 
data.  Within  the  inset,  X  =  FIEPP-ICE  data,  C  = 
ealeulated  isentrope. 

III.  ACCURACY  OF  ICE  TECHNIQUE 


Figure  4.  Ideal  eurrent  waveforms  for  shoek-free 
ICE  measurements  up  to  2.2  TPa.  The  sample 
thieknesses  are  2,  4,  6,  8  and  10  mm. 

HEPP-ICE  is  the  only  teehnique  eapable  of  providing 
sueh  long  eurrent  waveforms,  and  therefore  the  only 
teehnique  eapable  of  produeing  shoek-free  eonditions  at 
these  large  stresses  and  in  sueh  large  samples.  One 
obvious  advantage  of  these  large  sample  sizes  is  that  the 
relative  errors  in  thiekness  and  time  measurement  are 
redueed,  see  seetion  III. 

II,  HEPP-ICE  EOS  RESULTS 

The  LANE  HEPP-ICE  program  has  foeused  on 
developing  the  HEPP  teehniques  to  produee  high  stress 
EOS  data,  and  relatively  few  data  gathering  experiments 
have  been  performed.  The  data  for  OFHC  eopper  were 
presented  previously  [3].  Here  we  will  present  the  data 
for  pure  tungsten.  There  are  two  teehniques  we  use  for 
data  analysis,  the  generalized  Lagrangian  approaeh 
deseribed  in  the  Introduetion,  and  the  Baekward  teehnique 
[14].  Due  to  spaee  limitations,  only  the  results  of  the 
generalized  teehnique  are  shown. 

The  results  for  eopper  showed  that  the  HEPP-ICE  is 
eapable  of  resolving  differenees  between  isentropes  as 
small  as  0.2%.  These  EOS  data  for  pure  tungsten  were 
obtained  from  a  series  experiments  were  the  sample 
thieknesses  ranged  from  1  to  4  mm.  The  pressure  vs. 
partiele  veloeity  {P,  u)  isentropes  were  obtained  by  fitting 
quadraties  to  the  stress  vs.  u  data  above  the  elastie-plastie 
transition,  then  setting  the  eonstant  term  to  zero  so  that  the 
data  interseeted  the  origin;  see  Figure  5  whieh  shows  the 
data  and  the  ealeulated  isentrope.  The  data  agree  to 


The  aeeuraey  of  the  ICE  teehnique  is  limited  by  the 
aeeuraey  of  the  wave  and  partiele  veloeity  measurements, 
and  the  uniformity  of  the  stress  loading.  Substituting  A ’s 
for  the  J’s  in  Equation  (1)  to  represent  the  finite  steps  in 
the  data,  the  fraetional  error  of  steps  in  stress  is 


dAcr  dh  dAu  .  ,  dAt 

- » - 1 - he,  (m) - 

Act  h  Au  h 


(2) 


where  the  three  terms  on  the  right  hand  side  are  the 
fraetional  errors  in  thiekness,  partiele  veloeity  and  transit 
time.  Note  that  the  transit  time  is  At  =  h  ^  Ci(u),  so  the 
timing  error  inereases  with  wave  veloeity  Ci(u). 

The  inaeeuraey  of  the  VISAR  surfaee  veloeity 
measurement  ean  be  as  small  as  0.1%  with  eareful 
preparation,  and  may  be  ~2%  otherwise  [15].  The 
measurement  of  wave  veloeity  is  Cl(u)  =  h  ^  At,  and  both 
h  (thiekness)  and  At  (time  differenee)  are  sourees  of  error. 
Presently,  h  may  be  measured  to  an  aeeuraey  of  1  pm, 
whieh  represents  an  error  of  0.01%  for  a  10-mm  thiek 
sample. 

The  errors  in  time  are  eaused  by  the  laek  of 
synehronization  between  VISAR  signals  and  by  the  time 
resolution,  dt,  of  the  VISARs.  After  synehronization  of 
the  data  there  remains  a  possible  systematie  error  of  up  to 
+V2  dt  between  ehannels.  This  error  ean  have  a  signifieant 
effeet  on  the  aeeuraey  of  Cl-  For  example,  for  =  1  mm. 
Cl  =  7  km/s,  and  At  =  143.86  ns,  an  error  of  ±1/2  ns 
would  render  an  error  in  Cx  of  up  to  ±  0.35%. 

Despite  these  systematie  timing  errors,  for  the  results  to 
be  eonsistent  they  must  all  he  on  the  same  isentrope. 
Consequently,  the  data  may  be  adjusted  by  a  minimization 
algorithm  that  adds  a  fixed  time  eorreetion  to  eaeh  eurve 
in  turn  until  the  standard  deviation  between  the  Cl  vs. 
stress  eurves  is  minimized  and  the  data  are  self-eonsistent; 


1366 


Table  1.  Relative  errors  in  ICE  experiments. 


VI.  REFERENCES 


Error 

HEPP-ICE  error 

B-field  non-uniformity  aeross 
sample,  a  few  % 

<0.1% 

B-field  sample -to-sample  non¬ 
uniformity,  a  few  % 

Zero 

Thiekness:  ±0.1%  for  1-mm 

±0.01% 

VISAR  timing:  for  1-mm  at 
10km/s  =  ±0.1%,  ±0.5% 

±0.01%,  ±0.05% 

VISAR  veloeity  ±0.1%  to  ~2% 

Same 

these  adjustments  must  lie  within  ±'72  dt.  This  algorithm 
works  well  but  presently  laeks  mathematieal  rigor. 
However,  using  this  algorithm  the  standard  deviation 
between  reeords  was  redueed  to  ~0.2%  for  the  eopper 
data.  These  data  demonstrate  the  eapabilities  of  the 
teehnique,  but  until  we  ean  devise  a  rigorous  algorithm 
we  eannot  assert  that  0.2%  is  the  true  error.  The  use  of  a 
veloeimeter  (VISAR  or  PDV  [16])  with  better  time 
resolution  is  preferred.  All  errors  are  summarized  in 
Table  1. 

IV.  SUMMARY 

It  has  been  shown  that  the  HEPP-ICE  load 
eonfiguration  ean  produee  magnetie  stress  uniformity  to 
better  than  1  part  in  1000  over  the  eentral  87%  of  the 
sample  faees,  and  it  provides  exaet  matehing  of  the 
stresses  between  opposing  samples.  A  large-seale  HEPP- 
ICE  system  is  eapable  of  produeing  shoek-free  loading  up 
to  2.2  TPa  in  10  mm  thiek  tungsten  samples. 

The  aeeuraey  of  the  HEPP-ICE  system  is  enhaneed  by 
the  uniformity  and  matehing  of  magnetie  stress  loading, 
and  by  the  larger  sample  sizes  that  ean  used  in  HEPP-ICE. 
Data  for  tungsten  (and  eopper)  have  been  presented  whieh 
demonstrate  the  inaeeuraeies  of  the  teehnique,  whieh  may 
be  as  low  as  0.2%  in  stress. 
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